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CONSTRAINING PRIMORDIAL MAGNETIC FIELDS THROUGH LARGE SCALE STRUCTURE 
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ABSTRACT 

We study primordial magnetic field effects on the matter perturbations in the Universe. We do 
not limit analysis by considering a particular magnetogenesis scenario that assumes magnetic field 
generation prior to recombination. Contrary to previous studies, we limit the total magnetic field 
energy density and not the smoothed amplitude of the magnetic field at large (order of 1 Mpc) scales. 
We review several cosmological signatures, such as halos abundance, thermal Sunyaev Zel'dovich (tSZ) 
effect, and Lyman-a data. For a cross check we compare our limits with that obtained through the 
CMB faraday rotation effect and nucleosynthesis bound. The limits are ranging between 1.5 nG and 
4.5 nG for n B e (-3; -1.5). 

Subject headings: primordial magnetic fields; early universe; large scale structure 



1. INTRODUCTION 

Observations show that galaxies have magnetic fields 
with a component that is coherent over a large frac- 
tion of th e galaxy with field strength of order 10~ 6 
Gauss (G) (IBeck et al.lll^lWidrowll2fMlvTO[200l . 
These fields are supposed to be the result of amplifica- 
tion of initial weak seed fields of unknown nature. A 
recent study, based on the correlation of Faraday rota- 
tion measures and Mgll absorption lines (which trace 
halos of galaxies), indicates that coherent /iG-strength 
magnetic fields were already in place in normal galax- 
ies (like the Milky Way) when the universe w as less 
than half its present age ([Kronberg et al.l l2008f ) . This 
places strong constraints both on the strength of the 
initial magnetic seed field and the time-scale required 
for amplification. Understanding the origin and evolu- 
tion of these fields is one of the challenging questions of 
modern astrophysics. There are two generation scenarios 
under discussion currently: a bottom-up (astrophysical) 
one, where the needed seed field is generated on smaller 
scales; and, a top-down (cosmological) scenario, where 
the seed field is generated prior to galaxy formation in 
the early universe on scales that are large now. More 
precisely, astrophysical seed field sources include battery 
mechanisms, plasma processes, or simple transport of 
magnetic flux from compact systems (e.g. stars, AGNs), 
where magnetic field generat ion can be extremely fas t 
because of the rapid rotation ([Kulsrud fe Zwe ibcl 2008). 
Obviously, the correlation length of such a seed field 
cannot be larger than a characteristic galactic length 
scale, and is typically m uch smaller. In the cosmolog- 
ical seed field scenario, (|Kandus et al.l I2T5TTI ). the seed 
field correlation length could be significantly larger than 
the current Hubble radius, if it was generated by quan- 
tum fluctuations during inflation. There are different op- 
tions for seed field amplification, ranging from the MHD 
dynamo to the adiabatic compressio n of the magneti c 
field lines during structure formation ([Beck et al.lll996l ). 



The presence of turbulence in cosmic plasma plays a cru- 
cial role in both of these processes. The MHD turbu- 
lence was investigated a long time ago when consider- 
ing the processes in astrophysical plasma, while there 
is a lack of studies when addr essing the turb ulence ef- 
fects in cosmological contexts (Biskamp 2003j). At the 
late stages of the universe evolution the energy den- 
sity present in the form of turbulent motions in clus- 
ters can be as large as 5-10% of the thermal energy den- 
sity (jKravtsov fe Bprgani Il2012f). This can in fluence the 
physics of clusters ()Subramanian et al.ll2006l ). and/or at 
least should be mo deled correctly wh e n performing larg e 
scale simulations (|Vazza et al.l 120061 : iFeng et all [2009) . 
The proper accounting of the MHD t urbulence effects is 
still under discussions (jSprin gcl 2010). Both astrophysi- 
cal and primordial turbulence might have distinctive ob- 
servational signatures. As we already noted above, the 
most direct signature of MHD turbulence is the observed 
magnetic fields in clusters and galaxies. 

Galactic magnetic fields are usually me asured t hroug h 
the induced Faraday rotation effect (see iValled ([2004)) 
and, as mentioned above, the coherent field magnitude 
is of order a few fiG with a typical coherence scale of 
10 kpc. 1 On larger scales there have been recent claims 
of an observed lower limit of or der 10~ 15 — 10~ 16 G on 
the intergalactic magnetic field dNeronov fe Vovkl 120101 : 
iTavecchio et al.ll2010l : lDolag et al.ll2011| ). assuming a cor- 
relation length of A > 1 Mpc, or possib ly two orders 
of magnitude smaller ((Permer et all 1201 ID . An alterna- 
tive approach to exp lain the blazar s pectr a anomalies 
has been discussed bv lBroderick et ail (|2012f l . where two 
beam plasma instabilities were considered. 2 Although 
these instabilities are well tested through numerical ex- 

1 On the other hand, simulations starting from constant comov- 
ing magnetic fields of 10 _11 G show clusters generating fields suffi- 
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ciently large to explain Faraday rotat ion measurements Dol; 
p002ti ; IBaneriee iTJed amzik '( 2003 ) . 

2 The recent study Arlcn et al. ( 2012) claims that proper ac- 
counting for uncertainties of the source modeling leads to consis- 
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periments for laboratory plasma for a given set of pa- 
rameters such as a temperature and energy densities of 
beams and background, its efficiency might be questioned 
for cosmological plasma because of a significantly differ- 
ent (several orders of magnitudes) beam, and background 
temperature and energy densities. Prior to these ob- 
servations, the intergalactic magnetic field was limited 
only to be smaller than a few nG from cosmological ob- 
servations, such as the limits on the cosmic microwave 
back ground (CMB) radiat ion polarization plane rota- 
tion (jYamazaki ct al. 20f0) and on the Faraday rotation 
of polarized emis sion from distant blazars and quasars 
(jBlazi et al.lfl999h . 

In present paper we consider the presence of the pri- 
mordial magnetic field in the Universe and give a simpli- 
fied description of its effect on the large scale structure 
formation. We focus on the linear matter power spec- 
trum in order to show that even the total energy density 
presented in magnetic field (and as a consequence in mag- 
netized turbulence) is small enough, its effects might be 
substantial, and the effect becomes stronger due to non- 
linearity of processes under consideration. 

It has become conventional to derive the cosmologi- 
cal effects of a seed magnetic field by using its spec- 
tral shape (parameterized by the spectral index n B ) 
and the smoothed value of the magnetic field {B\) at 
a given scale A (wh i ch is usually taken to be I Mpc). In 
IKahniashvili et all (| 2 1 If) we developed a different and 
more adequate formalism based on the effective magnetic 
field value that is determined by the total energy density 
of the magnetic field. This formalism has been applied to 
describe two different effects of the primordial magnetic 
field; th e CMB Faraday rotatio n effect and mass dis- 
persion (|Kahniashvili et aDl2010f ). As a striking conse- 
quence, we show that even an extremely small smoothed 
magnetic field of 1 -29 G at I Mpc, with the Batchelor 
spectral shape (ng = 2) at large scales, can leave de- 
tectable signatures in CMB or LSS statistics. In present 
investigation we have been focus on the thermal Sunyaev- 
Zel'dovich effect, the cluster number densi ty, and Lyman- 
a dat a. These tests has been studied in Sh aw fc Lewisl 
(2010), but again in the context of the smoothed mag- 
netic field. We also do a more precise data analysis, and 
we focus on the inflation-generated magnetic fields. 

The structure of the paper is as follows. In Sec. II 
we briefly review the effective magnetic field formalism 
and discuss the effect on the density perturbations. In 
Sec. Ill we review observational consequences and derive 
the limits on primordial magnetic fields. Conclusions are 
given in Sec. IV. 

2. MODELING MAGNETIC FIELD INDUCED MATTER 
POWER SPECTRUM 

A stochastic Gaussian magnetic field is fully described 
by its two-point correlation function. For simplicity, we 
consider the case of the non-helical magnetic field, 3 for 
which th e two-point correlation fu nction in wavenumber 
space is (Kalmias hvili et aLll2010[) 

(Bt(k)B,(k')) = (27r) 3 S^(k~k')P ZJ (k)P B (k). (I) 

tence with a zero magnetic field hypothesis. 

3 We limit ourselves by considering the non-helical magnetic field 
because the density perturbations, and as a result the matter power 
spectrum, is not affected by the presence of magnetic hclicity. 



Here, i and j are spatial indices; i. j € (I, 2, 3), fe,- = ki/k 
is a unit wavevector; Pij(k) = Sij — kikj is the transverse 
plane projector; (S^ 3 ^(k — k') is the Dirac delta function, 
and P B {k) is the power spectrum of the magnetic field. 

The smoothed magnetic field B\ through the mean- 
square magnetic field, B\ 2 — (B(x) ■ B(x))|a, where the 
smoothing is done on a comoving length A with a Gaus- 
sian smoothing kernel function oc exp[— x 2 /X 2 ]. Corre- 
sponding to the smoothing length A is the smoothing 
wavenumber k\ = 2ir/\. The power spectrum P B (k) is 
assumed to depend on k as a simple power law function 
on large scales, k < ko (where k]j is the cutoff wavenum- 
ber), 



P B (k) = P BO k ni 



2v 2 \*B\ 
r(n B /2 + 3/2) 



(\ky 



(2) 



and assumed to vanish on small scales where k > kp. 

We define the effective magnetic field B c g through the 
magnetic energy density p B — B 2 s /(8ir). In terms of the 
smoothed field, the magnetic energy density is given by 



B 2 (k D \) n *+ 3 
87iT(n B /2 + 5/2) 



(3) 



and thus Bos = B x {k D \)( nB+ ^/ 2 / \jY(n B /2 + 5/2). For 
the scale-invariant spectrum n B = — 3 and i? c ff = B\ 
for all values of A. The scale-invariant case is the only 
case where the values of the effective and smoothed fields 
coincide. For the causal magnetic fields with n B = 2 
the smoothed magnetic field value is extremely small for 
moderate value of the magnetic field. 

We also need to determine the cut-off scale kp . We as- 
sume that the cut-off scale is determined by the Alfven 
wave damping scale ku ~ vaLs, where v a is the Alfven 
veloc i ty and Ls the Silk damping s cale (jJedamzik et al.l 
Il998t [Subramanian fc Barrowilf 9981 ). Such a description 
is more appropriate when we are dealing with a homoge- 
neous magnetic field, and the Alfven waves are the fluctu- 
ations Bi(x) with respect to a background homogeneous 
magnetic field B (|Bi| <C |B |). In the case of the 
stochastic magnet ic field we generali ze the Alfven veloc- 
ity definition from lMack et al.l (|20 02) . by referring to the 
analogy between the effective magnetic field and the ho- 
mogeneous magnetic field. Assuming that the Alfven ve- 
locity is determined by -B G ff, a simple computation gives 
the expression of ko in terms of B c r: 



IMpc" 



= 1.4, 



(2 7r )"B+3/ 1 



10" 7 G 



T(n B /2 + 5/2) V B eS 



(4) 



Here h is the Hubble constant in units of 100 km s 1 
Mpc -1 . 

Note that any primordial magnetic field should satisfy 
big bang nucleosythesis (BBN) bounds, i.e., the total en- 
ergy density of the magnetic field should not be greater 
than 10% of the radiation energy density. Assuming that 
the magnetic field energy density is not damped away 
by MHD processes, the BBN limit on the effective mag- 
netic field strength, B c s < 8.4 x I0 -7 G, while trans- 
ferred in terms of B\ the BBN bo unds result in extremely 
small values for causal fields, see (|Caprini fc Durreri[200l1: 
IKahniashvili et ahll2011|) . 
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Fig. 1. — The magnetic field matter power spectra for different 
values of B\ (a) and for different values of tij (b). 

The primordial magnetic field affects all three kinds 
of metric perturbations, scalar (density), vector (vortic- 
ity), and tensor (gravitational waves) modes through the 
Einstein equations. The primordial magnetic field gen- 
erates a matter perturbation power spectrum with a dif- 
ferent shape compared to the standard ACDM model. 
As we noted above in this paper we focus o n the mat- 
ter p e rturbations. As it ha s been shown by (Ki m et all 
119961 : iGopal fc Sethil 120051 ). the magnetic-field-induced 
matter power spectrum P(k) is oc k 4 for tib > —1.5 and 
oc k 2nB+T for riB < —1.5. This in turn affects the for- 
mation of rare objects like galaxy cluster s which sample 
the ex ponential tail of the mass function. Sh aw fc Lewis! 
(2010) studies in great detail the formation of the mag- 
netic field matter power spectrum through analytical de- 
scription, and provides the modified version of CAMB 
that includes the possibility of a non-zero magnetic field. 
We have used the CAMB code to determine the matter 
power spectra for a wide range of the magnetic field am- 
plitudes and spectral indices. These spectra are shown 
in Fig. ([T]). It is obvious that the matter power spectrum 
is sensitive to the values of the cosmological parameters: 
the Hubble constant in the units of 100 km/sec/Mpc, 
h, fijif, and Q&, as well as to the density parameter of 
each dark matter component, i.e., f2 c d m and Q, v (here, 
M, b, cdm, and v indices refer to matter, baryons, cold 
dark matter, and neutrinos respectively, and f2 is the 
density parameter. To generate the matter plot we as- 
sume the standard flat ACDM model with zero curva- 
ture, and we use the following cosmological parameters: 
Q b h 2 = 0.022, n c h 2 = 0.1125, and h = 0.71. For sim- 
plicity, we assume massless neutrinos with three genera- 
tions. 4 As we can see the increase of the smoothed field 

4 The standard ACDM model matter power spectrum 
Pacdm(^) assumes a close to scale-invariant (Harrison-Peebles- 
Yu-Zel'dovich) post-inflation energy density perturbation power 
spectrum Po(fc) oc k n , with n ~ 1. 
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Fig. 2. — The tSZ power spectrum predictions at 150 GHz vary- 
ing the primordial magnetic field model at fixed cosmological pa- 
rameters, most importantly trg = 0.8. These predi ctions are com- 
pared against t he recent upper limits from ACT IIDunklev et ail 
[201Cf) and SPT flReichardt et al.ll20Tlj ) at £ = 3000. The current 
upper limits on the tSZ amplitude at I = 3000 do not constrain 
the primordial magnetic field parameters B e g and no as well as 
other observations. 

amplitude results in the additional power spectrum shift 
to the left, while increasing the value of ng makes the 
vertical shift. As we can see the large-scale tail (small 
wavenumbers) of the matter power spectrum is unaf- 
fected by the presence of the magnetic field. Below we 
address some of effects induced by the presence of the 
magnetic field, especially on large scales. 

3. OBSERVATIONAL SIGNATURES 

Primordial magnetic fields can play a potentially im- 
portant role in the formation of the first large-scale struc- 
tures. 

3.1. The Thermal Sunyaev-Zel'dovich effect 

As demonstrated in lShaw fc Lewis! (|2010f ) the strength 
of the primordial magnetic field affects the growth of 
structure. The power spectrum of secondary anisotropics 
in the CMB caused by the thermal Sunyaev-Zel'dovich 
effect (tSZ) is a highly sensitive probe o f the growth of 
structure (e.g. IKomatsu fc Selj aid 12001 ). The tSZ an- 
gular power spectrum probes the distribution of galaxy 
clusters on the sky essentially out to any redshift. At 
I ~ 3000, half of the contribution to the SZ power 
spectrum comes from matter halos with masses greater 
than ~ 2 x 10 14 M W at redshifts less than z ~ 0.5, see 
IBattaglia et"afl (|2012h . 

All the previous work on how primordial magnetic 
fields affect the tSZ powe r spectrum used the model from 
IKomatsu fc Seliakl (|2002fl . here referred to as KS model, 
which has been shown t o be incompat i ble w ith recent 
observations of clusters (|Arnaud et al.l I2009D and tS Z 
power spectrum measurements by ILueker et al.l ([2009). 
Using the KS model for primordial magnetic field stud- 
ies also ignores all the recent advancements in tSZ power 
spectrum theory and predictions that illustrate the im- 
portance of properly modeling th e detailed astrophysics 
of th e intracluster med i um (e.g. IBattaglia et al.l 120101 
l2012t IShaw et all l2010t iTrac et all 120 111) . We modify 
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the code described in [Shaw fc Lewis! (|2010f ) to include 
these improvements by changing the pressure profile 
used in their model from KS to the profile given in 
iBattaglia et~aTl (|2010L l2012t ). The results from the new 
pressure profile are shown in Fig. ([2]) with the great- 
est difference being the amplitude of the new tSZ power 
spectrum is approximately two times lower than pre- 
vious predictions and below the current o bservations 
constraint from ACT (jDunklev et al.l l2010f ) and SPT 
(jReichardt et all 1201 If ) at I = 3000. Updating the the- 
ory predictions for the tSZ power significantly reduces 
the constraints put on primordial magnetic field param- 
eters using these observations. In Fig. @ we illustrate 
that the magnetic fields with an effective amplitude of 
order of 5 nG are almost excluded. Given that there 
is additional u ncertainty in the theoretical modeling 
of th e tSZ (e.g IBattaglia et al.l l2010L 120121 : iShaw et al.1 
2010; iTrac et alll2f)TTIV combined with si gnificant contri- 
butions from other secon dary sources (jReichardt et al.l 
[20111: iDunklev et al.ll20Toh around t ~ 3000, for example 
from dusty star forming galaxies, future tSZ power spec- 
trum measurements are not going to be competitive in 
constraining primordial magnetic fields parameters. 

3.2. Halos Number Density 

The predicted halo number density iV pre( j(M > Mq, z) 
depends on the considered cosmological model. One of 
important characteristics of a cosmological model is the 
linear matter power spectrum that we reviewed in Sec. 
II above. Below we discuss the halo number counts de- 
pendence on the presence of the magnetic field. 

The halo mass function at redshift z is N(M > 
M ,z) = J^ dM n(M,z), where n(M,z)dM is the co- 
moving number density of collapsed objects with mass 
lying in the interval (M,M + dM), and it can be ex- 
pressed as 

n(M,z) — -77^^/(^)377- (5) 



M 



dM 



The multiplicity functio n vf(v) is a universal func- 
tion of the peak height ([Press fc Schcchtcr 1974) v = 
5c/cr(R), where <j(R,z) the r.m.s. amplitude of den- 
sity fluctuations smoothed over a sphere of radius R — 
(3M/47rpj\/) 1 / 3 , and the critical density contrast 5q — 
1.686 is the density contrast for a linear overdensity able 
to collapse at the redshift z. Here, pm is the mean matter 
density. For gaussian flu ctuations vf{v) oc exp[— v 2 j2\ 
(|Press fe Schechter]|1974f ). where the normalization con- 
stant is fixed by the requirement t hat all of the mass 
lie in a given halo J vj{v)dv = 1/2 (White 2002). The 
evolution of the halo mass function n(M, z) is mostly 
determined by the z dependence of a{R, z). 

The r.m.s amplitude of density fluctuations <J 2 (R, z) 
is relate d to the linear matt er power spectrum P(k, z) 
through ([Jenkins et al.ll2001f) 



a 2 (R,z) 



D{zf 
2tt 2 



P{k,z)\W(kR)\ 2 k 2 dk, (6) 



where D(z) is the growth factor of linear perturba- 
tions normalized as D(z = 0) = 1 today, W{kR) is 
the Fourier transform of the top-hat window function, 
W{x) = 3(sina; — x cosx)/x 3 . In Fig. [3] we illustrate 
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Fig. 3. — v{M, z = 0) for different effective magnetic field values 
B e g and spectral index n-g. 

a(M, z — 0) function for the different values of the effec- 
tive magnetic field, B Q g, and the spectral index iib- The 
smaller amplitude of the magnetic field results in modi- 
fications at smaller mass scales. The a(M) dependence 
on the magnetic field ch aracteristics is also derived in 
([Kahniashvili et aDl2010f ). but contrary to the case pre- 
sented here, reflects only the a(M) induced by the pure 
magnetic field. In the present work we derive the effect 
from the magnetic field on the overall matter dispersion, 
including the standard density perturbations. The value 
of er(Af) at M = 2 x 10 14 M Sun is a round 0.8 that agrees 
well with observational data, see ([Burenin fc Vikhlininl 
[20T1 . 

Numerical computation results for n(M, z) are not ac- 
curately fit by the PS expression v f{v) cx exp[— ^ 2 /2] , 
see Refs. (ISheth fc Tolnrer] [l99l I Jenkins et all l200ll: 
IHu fc Kravtsovll2003f ). Several more accurate modifica- 
tions of n(M, z) have been proposed. H ere, we use the 
ST modificati onlSheth fc Tolmenl (|1999f ). as defined, (see 
Eq. 5 of Ref. (|Whitdl2002f) T~ 

f(v) cx [1 + (ais 2 )- p }(av 2 )- 1 / 2 exp[-av 2 /2} , (7) 

where the parameters a — 0.707, a nd p = 0.303 are fixed 
by fitting to the numeric al results (I White! 1200 ID (for the 
PS case: a = 1 andp = 0) ISheth fc Tolmenl ([1999T ). With 
this choice of parameter values the mass of collapsed ob- 
jects in Eq. ([7]) must be defined using a fixed over-density 
contrast with respect to the background density pu , and 
this requires accounting for the mass conversion between 
Afiso& and M200C Such a conve rsion depends on cosmo- 
logical parameters, see Fig. 1 of lWhitd l200l . Here, wo 
use an analytical extrapolation of this figure to do the 
conversion for Hm € (0.2,0.35). 

The difference induced by the magnetic field in the 
matter power spectrum P(k) can potentially modify 5c 
parameter entering in Eq. ([7|), that will result in differ- 
ent halo number counts. On the other hand, here we 
focus on the first order effects, so we neglect all changes 
induced by magnetic field in the Sheth-Tolm en model 
parameter fitting (see ISheth fc Tolme n (1999)). We also 
use the halo number count function at z — because we 
are focusing only on linear power spectrum, and all ef- 
fects relate^_Jo_Jhe_jna^etefield non-linear evolution 
(see [Schleicher & Miniay (|2011f) ) during the structure 
formation are neglected. We will present a more realistic 
scenario of the first object formation in future works. 

In Fig. ([4]) (top panel) we illustrate the halo mass func- 
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Fig. 4. — Halo number density N(M > Mo) (top panel) and 
ratio of number density for magnetic and non-magnetic simulations 
Nb)/No (bottom panel) for different effective magnetic field values 
-B e ff and spectral index «b, and 2 = 0. Number of small mass 
objects (M ~ 10 4 Mq) in magnetized case can be reduced down by 
factor of 100 compared to the non-magnetic number, object number 
count excess occurs for objects with mass around (M ~ 10 10 Mq). 

tion today (z = 0) for different values of and ub- 
As we can see, the magnetic field presence affects the 
small mass ranges, reducing the abundance of the low 
mass objects. We do not present here any statistics us- 
ing halo data accou nting for several uncer tainties involv- 
ing clusters physics (Battag lia et al.l l2012). On the other 
hand, we would like to underline that the presence of 
high enough magnetic field might be a possible explana- 
tion of the low mass objects abundance, which is one of 
the unsolved puzzles in ACDM cosmologies. 

To get a better understanding of the magnetic field in- 
fluence on the halo abundance, we plot the ratio of halo 
number density of ACDM models with and without mag- 
netic fields (see Fig. HI bottom panel). In the high mass 
limit all magnetized ACDM models predict slightly (a 
relative difference of the order of 10 -5 ) higher halo num- 
ber density. Number density excess peaks around halos 
with mass (M ~ IO 1O M0) and is strongly affected by 
the effective magnetic field value, as well as on the spec- 
tral shape. In contrast, at low mass limit M < 10 7 M Q , 
number of objects can be significantly lower as then its 
non-magnetic value. 

3.3. Lyman-a data 

The small scale modifications induced by the primor- 
dial magnetic field must be reflected in the first ob- 
ject formation in the Universe, i.e., the objects at high 
redshifts. The most important class of such objects 
are damped Lyman-a absorption systems. 5 To describe 

5 These objects have a high column density of neutral hydrogen 
(JV H i > 10 20 cm- 2 ) and are detected by means of absorption lines 
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Fig. 5. — The magnetic field m atter power s pectra for different 
values of ng and data points from Croft (2002) 

these systems it is possible to use semi-analytical mod- 
eling. Lyman-a systems has been used to constrain dif- 
feren t cosmological scenario, see Ref. ([McDonald et al.l 
2004), and references therein. Lyman-a data, which 
are very sensitive to the matter power spectrum around 
k ~ I0 -1 — IO 2 Mpc -1 , wavenum bers that are affecte d 
by the primordial magnetic field (|Shaw &: Lewis! [2010) . 
As we will see below these systems can be used to place 
stringent constraints on magnetic field properties. 

We do not go through the detailed modeling of Lyman- 
a systems, leaving this for more precise computations, 
but we use the direct comparison of the reconstructed 
matter power spectrum and the theoretical matter power 
spectra affected by the primordial magnetic field. 

For this study we use Ly man- a data obtained by the 
Keck telescope (|Croftll2002h . To get a conversion of data 
points (accounting that we use the wavevector k units 
/i/Mpc, we multiply data by the conversion factor 



100 



T+z 



given m Ref. (|Kim et alJl!996h . As the data is given at 
rcdshift 2.72, we translate the data to redshift zero by 
multiplying it by the square of the ratio of the growth 
factor at redshift zero to that at redshift 2.72. We com- 
pute the growth factors using the ICOSMOS calcula- 
tor. 6 Thus, we multiply the data by 8.145 to estimate 
the Lyman- data at redshift z = 0. The comparison 
of the theoretical predicted matter power spectrum and 
Lyman- data is given in Fig. ([5]). 

We use x 2 statistics to compare the predicted model 
with Lyman-a data. We assume no correlation between 
the uncertainties in the -P(k) measurements for different 
k values and find no evidence for primordial magnetic 
fields. 

The 95% and 68% confidence level limits are given in 



in quasar spectra (Wolfe 1993). Observations at high redshift have 
lead to estimates of the abundance of neu tral hydrogen in damped 
Lyman-a systems (Lanzctta et al. 1993). The standard view is 
that damped Lym an-a systems are a population of protogalactic 
disks (Wolfe 1993), with a minimum mass of M = 10 10 /i _1 M Slln 
(IHaehnerrJll995l ). 

6 ICOSMOS Calculator is available at 
http: / /www. icosmos.co.uk/index. html 
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Fig. 6. — The effective magnetic field limits from the lyaman-a 
data for different values of rig 
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Fig. 7. — The smoothed magnetic field limits from the Lyman-a 
for different values of ns 



Fig. |U The limits on B\ are given Fig. [7J We explic- 
itly present the limits for B c g and B\ just to show that 
they have different behaviors when the spectral index is 
increasing. In terms of the total energy density of the 
magnetic field the limits are weaker if we are considering 
the redder spectra. At this point the total energy density 
of the phase transition generated magnetic field is almost 
unconstrained. 

3.4. The CMB Faraday Rotation effect 

As we already noted above the primordial magnetic 
field induces the CMB polarization Faraday rotation, and 
for an homogen eous magnetic field the rotation angle is 
given through, (|Kosowskv & Loeblll996T ) 



a 



1.6° 



Bn 



InG 



30 GHz V 



(8) 



where Bq is the amplitude of the magnetic field, and v — 
is the frequency of the CMB photons. In the case of the 
stochastic magnetic field we have to determine the r.m.s. 
value of the rotation angle, a rms , and corresponding ex- 
pression in terms of the effective magnetic field is given 



Fig. 8. — The effective magnetic field values for different spectral 
index ng. Solid and dashed lines correspond to the 95% and 68% 
confidence levels, respectively. Upper limit set by the Big Bang 
nucleosynthesis £? e ff = 840 nG is shown by horizontal black line 
(BBN). 

in (jKahniashvili et al"]l2010[ ). being 



;0.14° 



D, 



err 



InG 



100 GHz 



V«s + 3 



y\(2i+i)i{i+i) / dx 

.1=0 Ja 



"n J (fc D 77o)("*+ 3 )/ 2 



x nB jf(x) 



(9) 



Here, rjo is today value of conformal time, ji(x) is a 
Bessel function with argument x — krjo, and xs — fcs^o 
where k$ = 2 Mpc -1 is the Silk damping scale. In the 
case of an extreme magnetic field which just satisfies the 
BBN bound, hp might become less than the Silk damp- 
ing scale. In this case the upper limit in the integral 
above must be replaced by xu = krji]Q. Note, that for 
iib — > — 3 Eq. (Q3) is reduced to Eq. ((8]) (see for details 
Ref. (jKahniashvili et al.ll2010D. 

Here, we quote Ref. (jKomatsu et al.ll2010f ) in order to 
determine the upper limits for the r.m.s. rotation angle. 
Adding the statistical and systematic errors in quadra- 
ture and averagi ng over WMAP, QUa D and BICEP (see 
for details Ref. ()Komatsu et alJfeOlOD ) with the inverse 
variance weighting, it was obtained a = —0.25° ± 0.58° 
at (68% CL), or -1.41° < a < 0.91° (95% CL). We 
obtain for the r.m.s. value (absolute) of the rotation an- 
gle |a rms .| < 0.477° and |a rms | < 0.997° (68% CL. and 
95% C.L., respectively) assuming gaussian statistics. In 
Fig. (J8]) we display the upper limits of the effective mag- 
netic fields using the rotation angle constraints quoted 
above. Note that these limits order of magnitude bet- 
ter th an obtained previously in Ref. (K ahniashvili et aLl 
2010) where we used the WMAP 7 years data alone. 
For almost scale-invariant magnetic field the limits are 
around 0.5 nG. As we can see for ns > —0.5 the BBN 
limits on the effective magnetic field strength are stronger 
than those coming from the CMB faraday rotation effect. 
The situation is completely different when determin- 
ing the limits for the smoothed magnetic field i?A=iMpc 
with rig > — 2 , which are extremely strong from BBN 
(jCaprini fc Durrerl [2fJ0lT : iKahniashvili et alj |2011[ ). and 
moderate in the case of the large scale structures or the 
CMB birefrigence, see above. 
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4. CONCLUSION 

In this paper we study the large-scale signatures of cos- 
mological magnetic fields, such as the thermal Sunyaev- 
Zel'dovich effect, halo number density, and Lyman-a 
data. Due to several uncertainties present in tSZ and 
halos abundance tests we find that Lyman-a measure- 
ments provide the tightest constraints on the primordial 
magnetic field energy density. We express these limits in 
terms of the effective value of the magnetic field, B c h . In 
the case of the scale invariant spectrum ub = — 3 these 
limits arc identical to limits on the smoothed magnetic 
field B\, (smoothed over a length scale A that is con- 
ventionally taken to be I Mpc). For a steep magnetic 
field with spectral index ns = 2 the difference between 
the limits derived in terms of the effective and smoothed 
field is several orders of magnitude. Also limits have dif- 
ferent behavior with incr easing of ns- At this po int, as 
we underlined previously (iKahniashv ili et aLll2010D using 
the smoothed magnetic field can result in some confusion: 
the smoothed magnetic field at I Mpc scales is extremely 
small, while the total energy density of the magnetic field 
is maximal allowed by BBN bounds. The small values 
of the magnetic fields for tib =2 (that corresponds to 
the phase transition generated magnetic fields) might be 
treate d as non-relevance on these fields. For example, in 
Ref. (|Shaw fc Lewis! I2010D it is claimed that the mag- 
netic field with the spectral in dex greater than -2.5 is 
exclu ded (IShaw fc Lewis 



Ref. (Kahnia shvili et al 



2010), while as it is shown in 



201 If ) the magnetic field with 



extremely small smoothed field value B\ at A = 1 Mpc 
order of 10 -29 Gauss with the spectral index ns = 2 
can leave observable traces on the CMB and large scale 
structure formation. The limits are ranging between 1.5 
nG and 4.5 nG for ns <E (—3;— 1.5). These limits arc 
comparable for those from the CMB polarization plane 
rotation. 

Not e when this paper was in final stage of preparation 
Ref. (jPandev fc Sethi! [2012D appeared showing that at 
the magnetic field can be strongly constrained by the 
first object formation, in particular through the Lyman- 
a data. 
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